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Abstract 

We obtain solutions to the coupled Schrodinger-Poisson equations. The solutions describe the 
evolution of cold dark matter density perturbations in an otherwise homogeneous expanding Fried¬ 
mann universe. We discuss the relationships between descriptions of cold dark matter in terms of 
a pressureless fluid, in terms of a wavefunction, of a classical scalar field, and a quantum scalar 
field. We identify the regimes where the various descriptions coincide and where they differ. 
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I. INTRODUCTION 


The identity of dark matter is among the most tantalizing questions in science today jl[. 
Fortunately we live in an era where a large number of observations directly or indirectly bear 
upon this question. Foremost among these are measurements of galactic rotation curves, 
observations of gravitational lensing by dark matter clumps on various scales, the cosmic 
microwave background anisotropy observations, and broad surveys of the visible matter 
distribution such as the Sloan Digital Sky Survey. 

The observations require the dark matter to be, in first approximation, cold and collision¬ 
less. Collisionless means that the main force acting on dark matter and therefore the main 
force by which dark matter manifests its presence is gravity. The dark matter may have 
in addition weak non-gravitational interactions but the observations are consistent with the 
absence of non-gravitational interactions. That the dark matter is cold means that its pri¬ 
mordial velocity dispersion is small. By primordial velocity dispersion we mean the velocity 
dispersion that the dark matter particles have even in the absence of density perturbations. 
An upper limit of order 10“®c on the primordial velocity dispersion follows from the require¬ 
ment that free streaming of the dark matter particles does not erase density perturbations 
on the smallest scales on which they are observed. 

Particle candidates for the dark matter are also meaningfully constrained by the require¬ 
ment that they must £t comfortably with the particles that are already known to exist, i.e. 
those described by the Standard Model. There are three broad categories of dark matter 
candidates that are thought to fit well into the existing scheme of particle physics: weakly 
interacting massive particles (WIMPs), axions or axion-like particles, and sterile neutrinos. 
WIMPs are motivated by supersymmetric extensions of the Standard Model. Their mass 
is typically 100 GeV, and their primordial velocity dispersion of order 10“^^c. Axions are 
motivated by the Peccei-Quinn solution of the strong CP problem, the puzzle within the 
Standard Model why the strong interactions are P and CP invariant. The axion mass is 
thought to be of order 10“® eV/c^ and the axion primordial velocity dispersion of order 
10“^^c. Sterile neutrinos have mass of order a few keV/c^ and primordial velocity disper¬ 
sion of order 10“®c, at the limit of what is allowed. For this reason, sterile neutrinos are 
sometimes called “warm dark matter”. Axions and WIMPs are definitely cold dark matter. 

Cold collisionless dark matter may be described in the linear regime of the evolution of 
density perturbations as a pressureless fluid. This is the description of cold dark matter 
in calculations of the cosmic microwave background anisotropies [^. Since cold dark mat¬ 
ter plays an important role in this context and the calculations agree very well with the 
observations, the pressureless fluid description has high credibility. To obtain the cosmic 
microwave background anisotropies a full general relativistic treatment is necessary be¬ 
cause the relevant evolution occurs in part on length scales of order the horizon. However, 
on length scales much less than the horizon (i.e. for wavevectors much larger than the Hub¬ 
ble rate) dark matter density perturbations are correctly described by Newtonian gravity. 
Linear Newtonian perturbation theory is simple, well understood and agrees with the 
general relativistic description on length scales much less than the horizon, where many of 
the interesting phenomena in large scale structure formation occur. It is therefore a very 
useful tool. 

L. Widrow and N. Kaiser pointed out that, on scales much less than the horizon, cold 
collisionless dark matter can be described by a wavefunction satisfying the Schrodinger- 
Poisson equations. As is discussed below, in Section HI, the wavefunction description is in 
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many ways more powerful than the pressureless fluid description. It allows the introduction 
of velocity dispersion whereas the fluid description allows none. It can be used to describe 
multi-streaming and caustics in the non-linear regime, whereas the fluid description breaks 
down in that regime. Indeed, Widrow and Kaiser carried out numerical simulations of 
structure formation using a wavefunction satisfying the Schrodinger equation, in lieu of N 
bodies satisfying Newton’s force law equation. Several such simulations have since been 
carried out [Gj. 

In the present paper, we reproduce the results of Newtonian linear perturbation theory us¬ 
ing a wavefunction solving the Schrodinger-Poisson (sometimes called Schrodinger-Newton) 
equations. As far as we are aware, this had not been done before although related work, also 
using the Schrodinger equation to analyze the growth of density perturbations in the early 
universe, can be found in Refs. One of our motivations is to show that the formalism 

does indeed work as expected. However, our main motivation is to prepare the ground for 
an in-depth study of the dynamical evolution of axion dark matter. 

The proposal that the dark matter may be axions originates with the papers of Ref. [l^ 
which showed that axions are copiously produced during the QCD phase transition, at a 
temperature of order 1 GeV. The estimate of the cosmological energy density of the axions 
thus produced was obtained by a simple classical treatment of the axion held. The axions are 
extremely weakly interacting and therefore collisionless. They are non-relativistic shortly 
after being produced and subsequently red-shifted by the expansion of the universe. Thus 


they are very cold, as was already mentioned. It was emphasized in Ref. [ll| that the axions 


produced during the QCD phase transition behave as cold dark matter on all scales much 
longer than the wavelength of the axion held, hereafter called the de Broglie wavelength. 
So, although axions are much lighter than WIMPs, they behave in many circumstances in 
the same way as WIMPs. However, there are diherences. 

Obviously, axions behave diherently from WIMPs on the length scale of their de Broglie 
wavelength. One manifestation of this is the existence of a Jeans’ length for axion dark 

1^ . The formula for the 


matter. This was originally pointed out by the authors of Refs. [12| 

Jeans’ length is given in Eq. fl2.3ip below. For QCD axions, those that solve the strong CP 
problem and have masses conservatively in the range 10“^ to 10“^^ eV, the Jeans length is 
far too short to affect structure formation in an observable way. However, we may also have 
axion-like particles (ALPs) with much smaller masses. Indeed, string theory predicts the 


15 . If the ALP mass is in the 10 
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existence of numerous axion and axion-like helds 
to 10“^"^ eV mass range and below, the Jeans’ length is large enough (kpc and larger) to 
have observable effects. Structure formation on length scales less than the Jeans’ length is 
suppressed. There is a long standing discrepancy by which observations show less structure 
on small scales than is predicted by N-body simulations. Many authors have proposed to 
resolve this by hypothesizing that the dark matter is an extremely light scalar held, with 
mass of order 10“^^ eV or less fl^ . 

Less obviously, axions differ from WIMPs because they thermalize and form a Bose- 
Einstein condensate (BEC) 0 0]. Axions thermalize as a result of their gravitational 
self-interactions when the photon temperature is of order 500 eV. Their thermalization time 
becomes shorter than the age of the universe then. When they thermalize, almost all axions 
go to the lowest energy state available to them. In this they differ from the other dark 
matter candidates. It was shown in Ref. T3| that, on all scales of observational interest. 


density perturbations in axion BEC behave in exactly the same way as those in ordinary 
cold dark matter provided the density perturbations are within the horizon and in the linear 


3 












regime. On the other hand, when density perturbations enter the horizon, and in second 
order of perturbation theory, axions generally behave differently from ordinary cold dark 
matter because the axions rethermalize so that the state most axions are in tracks the 
lowest energy available state. Axion BEG explains the occurrence of caustic rings of dark 
matter in galactic halos and their observed radii [l9[. It also solves the galactic angular 
momentum problem The observations require that at least 35% of the dark matter is 
axions j^ . 

Systems dominated by gravitational self-interactions are inherently unstable. In this 
regard the axion BEG differs from the BEGs that occur in superfluid ^He and dilute gases 
[2l|. The axion fluid is subject to the Jeans gravitational instability and this is so whether 
the axion fluid is a BEG or not jl^. The Jeans instability causes density perturbations to 
grow at a rate of order the Hubble rate H{t), i.e. on a time scale of order the age of the 
universe at the moment under consideration. Each mode of the axion fluid is Jeans unstable. 
However when the thermalization time is shorter than the age of the universe, the rate at 
which quanta of the axion held jump between modes is faster than the rate at which the 
Jeans instability develops. So the modes are essentially frozen on the time scale over which 
the axions thermalize. 

Our long term goal is to clarify the dual role of gravity in the evolution of the axion BEG. 
On the one hand gravity causes Jeans instability of the axion held modes. On the other, 
gravity causes axions to jump between those held modes. In the present paper, we take two 
steps towards this goal. In Section H we solve the Schrodinger-Poisson equations for self- 
gravitating collisionless dark matter. Our solutions describe the homogeneous expanding 
Friedmann universe and density perturbations therein. They also provide a complete set of 
states for the axions to occupy. In Section HI, we discuss the various relationships between 
descriptions of cold dark matter in terms of a pressureless huid, in terms of a wavefunction, of 
a classical scalar held, and a quantum scalar held, identifying the regimes where the various 
descriptions coincide and where they diher. In Section IV, we summarize our conclusions. 


II. WAVEFUNCTION DESCRIPTION OF LINEAR PERTURBATIONS 

Consider a huid composed of a huge number N of particles that are all in the same 
quantum-mechanical state. The wavefunction for the state satishes the Schrodinger 

equation 

idt'ip{f,t) = -h mT(r,t)j'0(r,t), (2.1) 

where m is the particle mass and T (r, t) is the Newtonian gravitational potential. In this 
section, we set h = c = 1. The density of particles in the huid is 

n(r, t) = V'0*(r, f)'0(r, t) . (2.2) 

Let us assume that the only kind of matter present is the huid of particles. The gravitational 
potential is then given by the Poisson equation 

= 47rGmn(r, f) . (2.3) 

The huid density satishes the continuity equation 

dtn + V -j = 0 , (2.4) 
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where 


(2.5) 


] = - i-^r) ■ 

The fluid velocity v{f,t) is deflned by = n{f,t)v{f,t). 

\/n{r, then 

v = —VP . 
m 

The velocity held satisfles the Euler-like equation 

dtv + {v ■ V)u = —VT — Vg , 


If we write = 

( 2 . 6 ) 

(2.7) 


where 


1 V^-v/n 
2m? y/n 


( 2 . 8 ) 


q is commonly referred to as “quantum pressure”. Eqs. (12.4^ and (12.7p follow from Eq. (I2.ip . 

We want to use Eqs. (12.ip and (12.dh to describe the evolution of density perturbations 
in an otherwise homogeneous Friedmann universe. The universe may be open or closed, 
or in between. However, because our description uses Newtonian gravity, the cosmological 
constant is set equal to zero. 

The wavefuction describing the homogeneous universe is 




(2.9) 


where H(t) is the Lemaitre-Hubble expansion rate. Indeed Eqs. (12.6p and (12.9p imply 


V = Hr 


( 2 . 10 ) 


The imaginary part of the Schrodinger equation is satisfied provided 


9*77,0 -h Shfuo = 0 , (2.11) 

and its real part is satisfied provided 

= + • ( 2 . 12 ) 

The Poisson equation then implies the acceleration equation 

+ =-^mno(t) . (2.13) 

3 

The continuity and acceleration equations may be combined as usual to yield the Friedmann 
equation 

, (2.14) 

a[t)^ 3 

where K = —1, 0, -|-1 depending on whether the universe is open, critical or closed, and a{t) 
is the scale factor deflned hy H{t) = ^. 

We now consider perturbations about this background: 


tfjif.t) = 7/10(7*, f) +7/11(7*, f). 


(2.15) 
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The perturbation is Fourier transformed in terms of comoving wavevector k as follows: 


= '0o(r,i) 


d^k "01 (fc, t)e'^ 


(2.16) 


Likewise the perturbation to the gravitational potential 




j d^k Ti(fc,f)e*^ . 


(2.17) 


The Schrodinger-Poisson equations expanded to linear order in the perturbations imply 


idt'ipi = + m(ToV'i + TiV’o), (2.18) 

2m 

and 

V'Ti = 47rGm(4>i + ^Po^Pl) • (2.19) 

It is useful to introduce the functions 

6{k,t) = ipi{k,t) + ipl{-k,t), (2.20) 

ri{k,t) = iJi{k,t)-'ipl{-k,t), (2.21) 

in terms of which Eqs. fl2.18p and fl2.19l) become 

idt6{k,t) - ^-^^7]{k,t) = 0, (2.22) 

(.. 3 , 

These can be combined into one, second order differential equation for 6{k,t): 

d^6{k, t) + 2H{t)dt5{k, t) - A7iGp5{k, t) + -^-^^6ik, t) = 0 , (2.24) 


where p = mriQ. The Fourier components of the perturbation to the wavefunction are given 
by 

t) = ^S{k, t) + dtS{k, t) (2.25) 

in terms of the solutions to Eq. (12.241) . 

The perturbation to the number density is 


ni{r,t) = 




4 k-r 

e , 


(2.26) 


and so the density contrast is 


6{r,t) 


no(r,t) 


J d^k 6{k, t)e^ 


(2.27) 
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By expanding 






one finds that 


Hence Eq. fl2.22p implies 


I3i{r,t) = — d^k r]{k,t)e-M 


^ X ia(t)k ^ . 

vi{k,t) = ^ ^ dt6{k,t) 
k ■ k 


(2.28) 

(2.29) 

(2.30) 


which is the same relationship between the velocity perturbation and the density contrast 
as in the standard description of cold dark matter in terms of a pressureless fluid. Eq. fl2.24p 
is also the standard second order differential equation governing the evolution of the density 
contrast, except for the last term. It arises due to quantum pressure in Eq. (12.7p and 
produces a Jeans length [H, 13 1 


^J = (lOvrGpm^) ^ = i.Ql ■ 10^'^cm 


10 °eV\2 ( 10 ^®g/cm' 


m 


P 


(2.31) 


For k > the Fourier components of the density perturbations oscillate in time. For 
k « the most general solution of Eq. fl2.24p is 

6(k, t) = A{k) + B(k) , (2.32) 

in the critical universe case [a{t) oc ts]. A{k) and B(k) are the amplitudes of growing and 
decaying modes, respectively. On distance scales much larger than the Jeans length, the 
wavefunction description coincides in all respects with the fluid description. 

Let us mention briefly that the wavefunction can also describe rotational modes, provided 
vortices are introduced. See, for example ref. [1^. In a region where V x E 7^ 0, the vortices 
have the direction of V x n and have density (number of vortices per unit area) ^|V x v\. 
By Kelvin’s theorem, the vortices must move with the fluid. Therefore, in an expanding 
universe, the density of vortices decreases as a(f)“^. Hence v oc a{t)~^ for rotational modes, 
which is again the usual result. 


III. DISCUSSION 

We saw in the previous section that density perturbations in the early universe may 
be described by a wavefunction which solves the Schrodinger-Poisson equations and that 
on length scales large compared to the Jeans length, Eq. fl2.3ip . the resulting description 
coincides with that in terms of a pressureless fluid. It is our purpose in the present section 
to place this result in a wider physical context. 

First let us state that, although it appears that the wavefunction description had not 
been explicitly given before, it is no surprise that it exists since the Schrodinger equation 
implies the continuity equation and the Euler-like equation fl2.7p . These two equations are 
the basic equations describing a fluid. The only difference is the quantum pressure term in 
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Eq. fl2.7p but that term is unimportant on distance scales large compared to the de Broglie 
wavelength. The Jeans length of Eq. fl2.3ip can be viewed as the de Broglie wavelength of 
the minimum energy state in a region of density p. Indeed in such a region, the gravitational 
potential is T = ^Gpr'^ and hence the energy of a trial wavefunction of width b is of order 

1 27r 

E(b) reaches its minimum for b ~ {^Gpm^)~^ ~ Ij. 

However the mathematical equivalence of the two descriptions hides important physical 
differences. This is perhaps best illustrated by an example. Consider the wavefunction 

. (3.2) 

where H is a constant and ^ It solves the Schrodinger equation for a free particle. 

The fluid with N particles in the state of wavefunction 'ip{x,t) has two flows, both with 
density ni = n 2 = iV|Hp, and with velocities vi = E and V 2 = — On the other hand, 
Eqs. fl2.2p and fl2.6p map onto a fluid whose density is n{r) = cos^(/c • r) and 

whose velocity n = 0. The two descriptions are mathematically equivalent in the sense that 
n{f,t) and v{f,t) satisfy Eqs. (12.21) and (12.7p because satishes Eq. (12.ip . But the two 

descriptions are not physically equivalent. They are physically equivalent only if we average 
over distances large compared to the wavelength ^ and if we ignore the velocity dispersion 
An = The spatial averaging is justihed in the limit k —)■ oo. Ignoring the velocity 
dispersion is justihed in the limit, ^ ^ 0. The two limits are compatible only if m —)• cx). 
This indicates that the physical differences between the wavefunction and huid description 
disappear completely only in the limit where the dark matter particle is very heavy. 

The huid description never allows velocity dispersion since the velocity held v{r, t) has a 
single value at every point. In contrast, the wavefunction description allows velocity disper¬ 
sion and multi-streaming. The wavefunction description is richer therefore. It can describe 
everything that a huid describes but the reverse is not true. Whether either description is 
correct depends on the situation at hand. 

Consider a dark matter particle with properties typical of a WIMP candidate: m ~ 100 
GeV, density today uq ~ 10“®/cm^, and primordial velocity dispersion today Svq ~ 10“^^. 
The de Broglie wavelength associated with the primordial velocity dispersion is of order 10“^ 
cm, much smaller than the average interparticle distance of order 5 m. The particles are 
highly non-degenerate therefore. Provided that their primordial velocity dispersion is in fact 
irrelevant to whatever phenomenon is under study (free streaming would be an exception 
since it is a direct result of primordial velocity dispersion), the particles can be described 
as a pressureless huid. They can also be described by a wavefunction. The wavefunction 
description will in almost all respects be equivalent to the pressureless huid description but, 
unlike the latter, it allows the inclusion of velocity dispersion and its associated ehects. 
The wavefunction description is also applicable to the non-linear regime, after shell crossing, 
when the huid description in terms of a single velocity held v breaks down. The wavefunction 
description is more powerful because it packs more information. The wavefunction varies on 
a length scale of order 10“^ cm in the example given. The huid description is far coarser. 

Next consider a dark matter candidate typical of axions or axion-like particles: spin 
zero, m ~ 10“® eV, density today uq ~ 10®/cm^, and primordial velocity dispersion today 














6vo ~ 10“^^. The de Broglie wavelength associated with the primordial velocity dispersion 
is of order 10^® cm. The axion fluid is highly degenerate. The average occupation number of 
those states that are occupied is huge, of order 10®^. This suggests that axion dark matter 
is well described by a classical scalar held ip{f,t). The remainder of this section considers 
whether this is so. 

A classical scalar held satishes the Klein-Gordon equation 

- + ujlip + = 0 (3.3) 

where is the covariant derivative of general relativity. We allow the presence of a self¬ 
interaction = —^}P^ in the action density. First, let us emphasize that the classical held 
theory, Eq. fl3.3p . has no notion of axion. The axion is the quantum of the quantized scalar 
held which we call <h(r, f). There is no more notion of axion in Eq. fl3.3p than there is a 
notion of photon in Maxwell’s equations. Also there is no notion of mass since the mass m 
is the energy of an axion, divided by c^. Henceforth, for the sake of clarity, we no longer set 
h and c equal to one. wq in Eq. (13.31) is not the axion mass but the oscillation frequency of 
small perturbations in the classical scalar held in the inhnite wavelength limit. 

In the Newtonian limit of general relativity, the metric is ~ 2\h, g^i = 0, 

9ij = Eq. fl3.3p becomes then: 

(02 _ ^ 2^2 ^ 2 )^ ^ A 3 _ ^^2 + VT ■ V + = 0 . (3.4) 

We obtain the non-relativistic limit of this equation by setting 

= \/2 Re[e“*‘^°V(G t)] (3.5) 


and neglecting T versus c^, dtcj) versus versus a;o4/, and dropping terms proportional 

to and which indeed oscillate so fast as to ehectively average to zero. Eq. fl3.4p 

becomes then 


idtCp- '^V+o l<^lV+ 2^*^ • 

ZOJq OCUo 

(3,6) 

To obtain the Schrodinger equation 


ihdtip = + V{f,t)ip , 

2m 

(3.7) 

substitute 


0(Gt) = \ 

V ^0 

(3.8) 

in Eq. (13.bp and set m = The potential energy in Eq. (13. 7p is given, 

mean held theory, by 

in the sense of 

(3.9) 


The Newtonian limit of Einstein’s equation is the Poisson equation, Eq. fl2.3p . The non¬ 
linear version of Schrodinger’s equation obtained by substituting Eq. fl3.9p into Eq. fl3.7p is 
commonly called the Gross-Pitaevskii equation. 

So the Schrodinger equation describes the dynamics of a classical scalar held in the non- 
relativistic limit. This result is not new of course. We reproduced it here to prepare the 
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ground for the actual question we want to discuss, namely whether dark matter axions (or 
axion-like particles) are described by the Schrodinger-Poisson equations. Clearly, if axions 
are described by a classical scalar held, the answer is yes as we have just seen. But the axion 
is a quantum held. It may behave like a classical held some of the time or perhaps even all 
the time, but this is something that has to be proved. It cannot be merely assumed. 

Inside a cubic box of volume V = with periodic boundary conditions, the quantum 
axion held may be expanded (see for example Ref. fl^ l 

+ at(f)e-sP-'"] , (3.10) 


where n = (ui, n 2 , ns) with Uk {k = 1, 2, 3) integers, pa = ^n, u = ji\/p-p + c^m?. The an 
and at are annihilation and creation operators satisfying canonical equal-time commutation 
relations: 

[an{t),a\'i'^)]= ^n,n' , [a^(t), (f)] = 0 . (3.11) 

The classical held limit is the limit where h —?• 0 with hM held hxed, where M is the 
quantum occupation number of the state described by a particular solution of the classical 
held equations. The Hamiltonian for the quantum held $(P, t) which satishes Eqs. fl3.6p and 


fl2.3p in the classical held limit is |18 


H = 


hw. 


n 


E 


n2 ) 


(3.12) 


where is the sum of two terms: 

ni,n 2 


n\^n2,n^^n4^ 


723,714 


4 7l3,7t4 _ A 713,714 A 

^ni,n2 ' ^9 ni,n2 


The hrst term 


A 


n3,n4, 


'-s ni,n 2 ‘^ni+n2,n3+n4 

is due to the A<h^ type self-interactions. The second term 

AnGrri^h , 




A n3,n4 

ni,n2 




''ni+n2,n3-|-n4 


+ 


|Pni-Pn3P |Pni-Pn4 


is due to the gravitational self-interactions. The Heisenberg equations of motion are 


1 1 - 
= — — [if, ^ ^ A-^^ 






(3.13) 

(3.14) 

(3.15) 

(3.16) 


712,723,774 


We may likewise expand the classical held 


2Uf;V 


(3.17) 


The Fourier components Afi{t) satisfy 


iia. = 


(3.18) 


7^2,723,7^4 
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Eqs. fl3.16p and flS.lSp look similar but, as we will see, their physical implications are different 
because the are operators whereas the are c-numbers. 

Let us dehne the operator 

, (3.19) 

i.e. the occupation number at time t of the state labeled n. It was shown in Ref. [l8| that 
the Hamiltonian of Eq. fl3.12p implies the following operator evolution equation 

Ml = 'y ^ 2 \MiMj{Mi + 1)(A4 + 1) ~ MiMk{Mi + ^){Mj + l)] 2'n'6[u}i-{-u}j —ujk —uji) . 

k,i,j=l 

(3.20) 

To remove unnecessary clutter, we replaced fij by j. The derivation of Eq. fl3.2Up assumes 
only that the energy dispersion 6u of the highly occupied states is much larger than the 
relaxation rate T = -. r is the relaxation time, dehned as the time scale over which the 
distribution {A/)} changes completely. When dui >> T, the system is said to be in the 
“particle kinetic” regime. The same derivation that yields Eq. fl3.20p but applied to the 
classical counterparts 

Nft{t) = At{t)Afi{t) (3.21) 

yields 

Ni= l\A^S\‘'[NiNjNi + NiNjNk-NiNkNi-NkNiN,]27T6{ui+Uj-Uk-ui) , (3.22) 

k,i,j=l 


again in the particle kinetic regime. Eqs. fl3.20p and fl3.22p are clearly different, and they 
imply different outcomes. 

Consider the process i + j ^ k + I where two quanta, initially in states i and j move to 
states k and 1. Assuming ^ 0, this process always occurs in the quantum theory when 
the initial states are occupied. In the classical theory, the corresponding process occurs only 
if, in addition, at least one of the hnal states is occupied. In particular the scattering of 
two waves does not happen in the classical theory of Eq. fl3.6p if the only waves present are 
the two waves in the initial state. The quantum theory behaves differently because the hnal 
state oscillators have zero point oscillations. Incidentally, this observation shows that the 
oft repeated statement that the quantum and classical theories differ only by loop effects is 
incorrect. 

After a sufficiently long time, the classical and quantum systems thermalize and reach 
an equilibrium distribution. The time scale of thermalization of the classical system is 
of the same order of magnitude as that of the quantum system [l^ but the outcomes of 
thermalization are different. In the quantum case, Eq. fl3.2Up implies that the equilibrium 
distribution {Mj} is such that 

MiMj{Mi + 1)(A4 + 1) ~ MiMk{Mi + l)(A/)' + 1) = 0 (3.23) 


for every quartet of states such that Ui + Uj = Uk + ooi- Let us call 
Eq. (Km as 


whenever et + ej 


ek + ei. Eq. is solved by 



e 


huj, and rewrite 


(3.24) 


= C\n 



(3.25) 
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where C is a constant. Upon identifying C = this is seen to be the Bose-Einstein 

distribntion 


M = - . 

(3.26) 

QkgT _ 1 

On the other hand Eq. (|3.22p implies that the equilibrium distribution |Ah| for the classical 

case satishes 

(W + N,)NkNi = NiNj{Nk + Ni) 

(3.27) 

whenever e,- -1- Cj = e^. -1- ei. Eq. (13.271). which mav be rewritten as 

fill 

(3.28) 

is solved by 

M = c- . 

(3.29) 

Upon identifying C = ksT, we have 

Wu = ksT , 

(3.30) 


which is indeed the standard resnlt for classical oscillators at temperatnre T: each oscillator 
has energy ksT on average. 

We conclude that axion dark matter is not described by a classical held when it ther- 
malizes. Interactions are seen to have a dual role. They determine the behaviour of the 
classical held as described by Eq. fl3.6p . or Eqs. fl3.18p and fl3.22p which follow directly from 
Eq. 03.61) . Eq. 03. 6 p has a set of solutions which we may label The solutions de¬ 

scribe the states that the axions may occupy in the quantum theory. In the quantum theory, 
however, the interactions have the additional role of causing transitions between the various 
states Indeed if there were no such transitions the outcomes of thermalization in 

the quantum and classical theories would be the same. We just saw that they are not. 

Let r be the time scale over which the distribution {A/^} of the axions, over the states 
described by the classical solutions ^^(r, f), changes completely. We call r the relaxation or 
thermalization time scale. Note that full thermalization only happens generally on a time 
scale much longer than r. (The time scale for “full” thermalization depends on the degree 
of thermalization required and is therefore less robustly dehned than r.) On time scales 
short compared to r, axion dark matter behaves as a classical held because only relatively 
few transitions take place between the states described by the classical solutions ^^(r, t). 
On time scales long compared to r, the axions are not described by a classical held because 
their distribution over those states changes completely. On time scales long compared to r 
dark matter axions form a Bose-Einstein condensate (BEG) since they are highly degenerate 
and their number is conserved. The time scale for BEG formation is the relaxation time r 
18l. l23M25j . Axion BEG means that almost all axions go the lowest energy state available. 


The question is then: does r ever become shorter than the age of the universe t at that 
moment? It was found in Ref. fl^ . [l^ that T t during the QGD phase transition at a 
temperature of order 1 GeV when cold dark matter axions are hrst produced. The axions 
thermalize briehy then as a result of their A<I>^ interactions. Here and elsewhere, by the word 
‘thermalize’, we mean that the axion distribution relaxes and begins to approach a thermal 
distribution. We do not mean that they thermalize fully. This brief period of thermalization 
has no known observational consequences. However, when the photon temperature reaches 
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of order 500 eV, cold dark matter axions thermalize anew as a result of their gravitational 
self-interactions and this does have observational consequences, as was already mentioned 
in the Introduction. 


IV. SUMMARY 


We derived solutions of the coupled Schrodinger and Poisson equations. The solutions 
describe the homogeneous expanding matter-dominated universe and density perturbations 
therein. The description is identical to that obtained by treating the dark matter as a pres¬ 
sureless fluid, on all scales much larger than the de Broglie wavelength of the wavefunction. 
In a number of respects, the wavefunction description is simpler and hence superior. It has 
fewer degrees of freedom since a wavefunction is two real fields whereas a fluid is described 
by four real helds, the density and the three components of velocity. The mathematics is 
simpler as well. Even though it has only two real helds, the wavefunction can describe 
rotational modes. On the other hand, the meaning of the wavefunction is less intuitively 
obvious. 

We considered whether the wavefunction and huid descriptions are equivalent in general. 
They are equivalent in the sense that the density and velocity helds, given by Eqs. fl2.2p 
and fl2.6p . satisfy the huid equations if the wavefunction satishes the Schrodinger equation. 
However, the wavefunction and the huid describe objects which in general are not physically 
the same. In particular the wavefunction description allows velocity dispersion and multi¬ 
streaming whereas the huid description does not. The physical distinction between the two 
descriptions disappears completely only in the limit where the mass of the dark matter 
particle goes to inhnity. Because WIMPs are relatively heavy, the wavefunction and huid 
descriptions of WIMP dark matter are equivalent whenever velocity dispersion does not play 
a role. The wavefunction description has the advantage that it can describe phenomena 
associated with velocity dispersion, such as free-streaming, and that it can be used not only 
in the linear regime but also in the non-linear regime, after shell crossing. 

We asked whether the Schrodinger-Poisson equations correctly describe axion dark mat¬ 
ter. The answer is yes on time scales short compared to the relaxation time scale r, and 
no on time scales long compared to r. Whenever r is shorter than the age of the uni¬ 
verse t, axion dark matter is not correctly described by the Poisson-Schrodinger equations. 
Indeed axions move towards a Bose-Einstein distribution on the time scale r whereas the 
Schrodinger-Poisson equations would predict that they move towards a Boltzmann distribu¬ 
tion. Interactions, such as gravity or A<I>^ interactions, play a dual role. On the one hand 
the interaction influences the evolution of the classical held. The solutions of the classi¬ 
cal held equation, which is equivalent to the Schrodinger equation in the non-relativistic 
limit, describe the quantum states that the axions may occupy. But the interaction has the 
additional role of causing the axions to jump between those states. On time scales large 
compared to r, the distribution of quanta over the states described by the classical held 


changes completely. It was shown in Refs. j 171. IlSj that, as a result of axion gravitational 


self-interactions, r becomes and remains shorter than the age of the universe after the photon 
temperature has reached approximately 500 eV. Therefore the commonly made assumption 
that axion dark matter is adequately described by classical held equations at all times is 
incorrect. 

In summary then, we have obtained the solutions of the Schrodinger-Poisson equations 
that describe the homogeneous expanding Friedmann matter-dominated universe and den- 
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sity perturbations therein. Each solution corresponds to a possible quantum-mechanical 
state of dark matter axions. However the gravitational and other self-interactions of the 
axions cause them to jump between those states. The resulting dynamical evolution and 
observational implications of axion dark matter is the object of future work j^ . 


Acknowledgments 

P.S. would like to thank Sankha Chakrabarty and Yaqi Han for useful discussions. This 
work was supported in part by the U.S. Department of Energy under grant DE-FG02- 
97ER41209 at the University of Florida. 


14 


[1] G. Bertone, D. Hooper and J. Silk, Phys. Rep. 405 (2005) 279. 

[2] S. Dodelson, Modern Cosmology, Academic Press, Elsevier 2003. 

[3] E.M. Lifshitz, J. Phys. (Moscow) 10 (1946) 116; E.M. Lifshitz and I.M. Kalatnikov, Adv. 
Phys. 12 (1963) 185; J.M. Bardeen, Phys. Rev. D22 (1980) 1882; J.M. Bardeen, P.J. Stein- 
hardt and M. Turner, Phys. Rev. D28 (1983) 679; V.F. Mukhanov. H.A. Feldman and R.H. 
Brandenberger, Phys. Rep. 215 (1992) 203; K.A. Malik and D. Wands, Phys. Rep. 475 (2009) 

1 . 

[4] P. J.E. Peebles, The Large-Scale Structure of the Universe, Princeton U Press, 1980. 

[5] L.M. Widrow and N. Kaiser, Astrophys. J. 416 (1993) L71. 

[6] G. Davies and L.M. Widrow, Astroph. J. 485 (1997) 484; P. Coles and K. Spencer, M.N.R.A.S. 
342 (2003) 176; C.J. Short and P. Coles, JCAP 0612 (2006) 012, and JCAP 0612 (2006) 016. 

[7] I.M. Moroz, R. Penrose and P. Tod, Class. Quant. Grav. 15 (1998) 2733. 

[8] I. Szapudi and N. Kaiser, Ap.J. 583 (2003) LI. 

[9] R. Johnston, A.N. Lasenby and M.P. Hobson, MNRAS 402 (2010) 2491. 

[10] J. Preskill, M. Wise and F. Wilczek, Phys. Lett. B120 (1983) 127; L. Abbott and P. Sikivie, 
Phys. Lett. B120 (1983) 133; M. Dine and W. Fischler, Phys. Lett. B120 (1983) 137. 

[11] J. Ipser and P. Sikivie, Phys. Rev. Lett. 50 (1983) 925. 

[12] M.Y. Khlopov, B.A. Malomed and Y.B. Zeldovich, MNRAS 215 (1985) 575. 

[13] M. Bianchi, D. Grasso and R. Ruffini, Astron. Astrophys. 231 (1990) 301. 

[14] P. Svrcek and E. Witten , JHEP 0606 (2006) 051. 

[15] A. Arvanitaki et ah, Phys. Rev. D81 (2010) 123530. 

[16] S.-J. Sin, Phys. Rev. D50 (1994) 3650; J. Goodman, New Astronomy Reviews 5 (2000) 103; 
W. Hu, R. Barkana and A. Gruzinov, Phys. Rev. Lett. 85 (2000) 1158; E.W. Mielke and J.A. 
Velez Perez, Phys. Lett. B671 (2009) 174; J.-W. Lee and S. Lim, JCAP 1001 (2010) 007; A. 
Lundgren, M. Bondarescu, R. Bondarescu and J. Balakrishna, Ap. J. 715 (2010) L35; D.J. 
Marsh and P.G. Ferreira, Phys. Rev. D82 (2010) 103528; V. Lora et ah, JCAP 02 (2012) 011; 
T. Rindler-Daller and P. Shapiro, MNRAS 422 (2012) 135; D.J. Marsh and J. Silk, MNRAS 
437 (2014) 2652; B. Bozek, D.J. Marsh, J. Silk and R. Wyse, arXiv:1409.3544. 

[17] P. Sikivie and Q. Yang, Phys. Rev. Lett. 103 (2009) 111301. 

[18] O. Erken, P. Sikivie, H. Tam and Q. Yang, Phys. Rev. D85 (2012) 063520. 

[19] P. Sikivie, Phys. Lett. B 695 (2011) 22. 

[20] N. Banik and P. Sikivie, Phys. Rev. D88 (2013) 123517. 

[21] C.J. Pethik and H. Smith, Bose-Einstein Condensation in Dilute Gases, Cambridge University 
Press 2002. 

[22] N. Banik and P. Sikivie, arXiv:1501.05913. 

[23] D.V. Semikoz and LI. Tkachev, Phys. Rev. D55 (1997) 489. 

[24] K. Saikawa and M. Yamaguchi, Phys. Rev. D87 (2013) 085010. 

[25] J. Berges and J. Jaeckel, Phys. Rev. D91 (2015) 025020. 

[26] N. Banik, A. Christopherson, P. Sikivie and E. Todarello, in preparation. 


15 



